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Introduction 

A  complete  conception  of  large-scale  atmospheric  process-: 
determining  the  nature  of  the  weather  in  many  respects,  can 
only  be  obtained  by  the  inclusion  of  all  basic  factors  in  the 
baroclynic  atmosphere. 

In  the  analysis  of  large-scale  atmospheric  processes  we 
find  it  necessary  to  formulate  mathematical  theories  which  s) 
the  nature  and  evolution  of  the  real  fields  of  meteorological 
elements  as  exactly  as  possible.  In  principle,  the  mathematic 
constructions  based  on  the  application  of  the  laws  of  hydro- 
thermodynamics  in  atmospheric  conditions  permit  the  pertinent 
system  of  nonlinear  equations  to  determine  the  meteorologies,! 
elements  of  interest  to  us.  Therefore,  in  forecasting  meteor; 
logical  elements  the  solution  of  the  complete  system  of  honld 
hear  hydrothermodynamic  equations  appears  to  be  the ‘fundameir 
assignment  for  meteorological  theory. 


The  attempt  to  obtain  weather  forecasts  the  theoretical  way 
by  studying  the  system  of  hydrodynamic  equations  was  made  for 
the  first  time  by  Richardson  (9)  in  1922.  Richardson  did  r-cu,  c  ■ 
tain  sufficient  results  because  he  did  not  tackle  the  problem 
the  right  methodical  way.  Richardson  took  the  initial  wind  fiel< 
as  the  starting  feet  and  hoped  to  obtain  the  evolution  of  at¬ 
mospheric  currents  with  time  on  the  basis  of  the  solution  of  th 
equation  system.  Inasmuch  as  the  initial  field  of  velocity  gar.v 
rally  includes  manifold  perturbations  of  all  kinds  of  scales, 
Richardson  should  have  taken  into  account  the  evolution  of  all 
types  of  small  (sound,  gravitational)  perturbations  together 
.  with  the  evolution  of  large-scale  perturbations  which  are  of 
meteorological  interest,  although  the  former  are  insignificant 
in  meteorological  respects.  The  description  of  the  small  per¬ 
turbations  require  the  utilisation  of  correspondingly  small 
scales  and  small  time  integrals.  The  length  and  time  scales  us 
by  Richardson,  were  much  larger  than  the  degree  of  exactitude 
required;  therefore,  the  sole  forecast  calculated  by  Rxchardso. 
did  not  stand  the  test. 

The  first  feasible  solution  of  the  equation  system  for  at¬ 
mospheric  processes  was  found  by  I* A,  Kibel  (2)  in  1940.  He 
S  based  his  analysis  on  the  assumption  that  large-scale  atmos¬ 
pheric  processes  present  themselves  as  quasi-geostrophic  pro¬ 
cesses*  This  permitted  I. A.  Kibel  to  take  the  initial  fields 
of  oressure  and  temperature  as  starting  factors  and  not  the 
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initial  wind  field;  his  problem-' was  the  calculation  of  changing 
pressure  and  temperature  fields  in  time. 

.  ••  However,  in  I. A.  Kibel’s  analysis  the  possibilities '  of  o;u," 
geostr.ophic  samples  were,  not  completely  exhausted.  Therefore, 
the  temperature  derivative  of  time  ■*«■«  was  only  determined  m 

c?  u 

a  calculation  of  the  horizontal  mass  transfer,  and  in  the  pres-' 

sure  change  the  convergence  factor  ^nd  the  divergence  of 
a  t 

the  air  in  the '-atmosphere  were  not  considered. 

Is  other  meteorological  analyses  of  the  Friedmann-Kotshin, 
School,  the  system  of  hybr ©thermodynamic  equations  was  used  mo:, 
efficiently;  the  vortex  velocity ' equation  was  examined,  as  a  pr 
-  mot ing'  equation  connecting  the  complete  alteration  of  the  vor¬ 
tex  with  the  divergence  of  the  air  on  the  horizontal  plane. • 

In  1941  M.B.  Shvetz  (6)  obtained  tbs  expression  for  the  ve 
tical  velocity  in  the  shape  of  the  integral  according  to  alti¬ 
tude  from  the  individual  derivative  of  the  vortex  velocity. 

In  1943  E.&.  Blinova  (1)  managed  to  obtain  the  prognostic 
equations  for  pressure:  and  temperature  T  on  a  certain  "average’ 
atmospheric  level  through  linearization, taking  into  account  the 
sphericity  of  the  surface'  of  the  earth.  E.N.  Blinova  started 


with  the  condition  of  maintenance  of  vorticity  of 


M  /Lj. 


on  the  mean  level  and  with  the  assumption  that  tempera¬ 
ture  changes  in  the  atmosphere  depend  only  upon  the  horizontal; 


transfer  of  homogeneous  masses  of  different  temperatures.  E.N* 
Blinova  integrated  the  system  of  linear  equations  obtained  in 
this  manner,  under  general  conditions  of  the  elementary  field- 
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Foreign  meteorologists  see  the  possibility  of  utilising  tiy- 
cr ©dynamic  equations  for  weather  forecasts  in  a  different  way, 
H,  Ertel  (8)  analyzed  the  equation  1 

>/.,P  1  ,  . 

H  Lf 

obtained  by  him  starting  with  the  nature  of  the  area  of  depen¬ 
dence  of  eouation  upon  the  field  of  function  (p,^p)»  bit; 

•  *  $  t 

result  was  a  meteorologically  incorrect  conclusion  according 
to  which  it  is  impossible  to  calculate  the  future  pressure  raa-' 
in  advance .  The  deficiency  was  clearly  reflected  in  the  result; 
of  A.M.  Obuchov,d  (3)  analysis  published,  in  1949  in  which  the 
author  obtained  a  more  complete  equation  in  the  examination  o i 
the  barotrope  sample  of  the  atmosphere: 
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H0  -  altitude  of  homogeneous  atmosphere;  the  so 
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We  will  keep  to  the  symbols  used  in  the  first  part  of  this 

•n  vv  “t  f'c  ’•}  G 

C.j-  t  y  vj  X  y,f  9  „ 


J  2  2 

v  x  —  y 

r  s  .1—  ,  K  (r)  -  Macdonald  function  of  sero  order. 

T  o 

In  1950  N.I*  Bulajev  obtained  ^ 

for  the  change  of  pressure  in  the  baroclynic  atmos- 

*tke  (.crt<>itiwJ  4  dcC£C*1$£  wttA  h6.iffkt^ 

phare  under/tlinear  temperature^  The  equation  for  the  pressure 
change  on  the  mean  level  of  the  troposphere  was  obtained  as 
follows: 


Or  Or  ,  .  .  «££ x 

— -  -  ki  -  k0(p,Ap)  ~  K-  (i  f  A  *.) 

1  at  2  ■  3 


T  »  average  temperature  of  troposphere;  -  a  cert- 
constant' quantity. 

This  analysis  also  brought  us  the  formula  for  the  calcu¬ 
lation  of  vertical  velocity  at  different  levels  of  the  atmor 
phare : 


ObJ  g 


a  (s)AfT%  p)  -  a  (a)  {(T,  dp)  f  (p,AT)J 


p  -  pressure  at  the  mean  level  of  atmosphere;  CL  C  s) ,  Oi  (z)  » 

■JL 

the  coefficient  dependent  on  the  altitude, 

Gharney * s  (?)  analysis  examined  the  barotropic  sample  of 
the  atmosphere  and  the  simple  baroclynic  sample.  The  author 
uses  the  equation  for  the  winds  in  geostrophic  approximatior 
and  the  equation  of  the  heat  flux  without  considering  the  Vi 
tical  current.  Gharney  obtained  Po&asson's ' equation  for  the 
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derivative  of  pressure  at  different  levels.  The  resulting 
equation  in  the  barctropic  case  was  integrated  according  to 
the  method  of  finite  difference  with  the  aid  of  electronic 
equipment.  To  integrate  the  equation,  Charney  was  forced  t 
supplement  artificial  boundary  conditions,  pressure  values 

/  ,-1  A-'  c  f 

and  Laplaoe^W/ipressure  as  time  functions  at  the  boundary  of 
the  area. 

The  baroclynic  atmosphere  sample  was  examined  in  the  ana¬ 
lysis  published  by  I. A.  Kibe.1  (1953)*  The  author  started  wit 
the  equation  of  the  vortex  in  geostropbical  approximation  a; 
with  the  equation  of  the  heat  influx  without  considering  the 
vertical  current,  I. A,  Kibel  integrated  Pouasson’s  equation 

for  the  desired  function  on  surfaces  p  ~  const  and  ob~ 

t 

Paired  the  solution  as  follov/s: 

2v;  r. 


1>  z 

*•>  df  V* 


1 

2f 


r,  v  u 

F(r,  y'  ,  p)  In  -=-  rdrd  f  - - , 

r  <lt 


'd\ 


F 


f  ,  p}  -  the  known  function  of  the  fields  of  meteorology 


?  z 


cal  elements  in  the  entire  atmosnhere.  -  the  mean  value 

'  ^ 

of  the  desired  value  at  level  p  up  to  the  circumference  of 

radius  rn .  Value  with  sufficiently  large  r  (1,000km! 

1  d w  1 

appears  small  and  may  be  eliminated.  Excluding  the  average 

value  of  function  F ( r , ^  ,  p)  under  the  integral  symbol  and 

r, 

integrating  the  function  of  influence  I. A,  Kibel  ob~  * 

r 

Pained  the  prognostic  formula 
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r. 


w  ~  *— *  «  For  purposes  of  simplification  it  was  assumed  that 
function  equals  F{r?  sp5  in  regard  to  the  examined 

_pO  .Liltt*  $ 

In  this  analysis  the  authors  intended  to  obtain  general 
.equations  for  pressure  and  temperature  -  changes  and  an  equa¬ 
tion  for  the  vertical 'currents-  which  takes  all  basic  factors 
of  the  baroclyaic  atmosphere  into  account.  Compared'  to'  pre¬ 
vious  analyses  the  complete  inclusion  of  dynamic  factors  in 
the  temperature  changes  appears  to  be  a  new  factor. 

In  the  integration  of  the  obtained  differential  equations 

the  authors  found  solutions  which  were  expressed  by  space 

,  vktti  ij 

integrals  mi  determinate* expression^ /»dependfeig  upon  the  dis¬ 
tribution  of  pressure  fields  and  temperature  in  space  as 
well  as  tha  congruent  influence  functions. 


1. Formulation  of  Problem 


arjwwi. 


.In  the  study  of  the  dynamics 'of  pressure  and  temperature 
changes  and  the  formation  of -vertical  motion  in  the  baroclyn 
atmosphere  we  used  the  system  of  hydrothenaodynamic  equation 
in  regard  to  -  atmospheric  processes  on  a  large  scale.  A  chare. 


teristic  feature  of  this  system  which  distinguishes  it  from 
conventional  hydrodynamic  equations,,  appears  to  be  the  exi¬ 
stence  of  the  deflective  force  of  the  rotation  of  the  earth 
in  the  equations  of  motion. 

The  local  study  of  space  regions  with  horizontal  'dimen¬ 
sions  over  several  thousand  kilometers  permits  us  to  consider 
the  surface  of  the  earth  as  a  level  surface  subject  to  the 
limitations  of  the  region  under  study  and  the  use  of  the  rec 
t angular  coordinate  system  in  the  initial  equations* 

The  analysis  excludes  regions  in  the  immediate  proximity 
of  the  equator  ( -  0-30°)  because  the  nature  of  the  motion 
shows  fundamental  changes  in  these  regions. 

The  inital  system  of  hydrodynamic  equations  is  assumed 
as  follows: 

Equations  of  Motion 
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Equation  of  Statistics 


<5  P 
da' 


(1 


(3 


Equation  of  Continuity 


dt  gf  d  tr  cp 


Equation  of  State 

P  s  f  RT.  0 

The  following  symbols  were  used  in  the  equations  (1.1)  \ 
(1.65: 

u,  v,vi  ~  components  of  velocity  vector  V  to  coordinate  axis 
p  «  pressure,  p  ~  density,  T  «  air  temperature,  2 coco,. 
to-  angular  velocity  of  earth  rotation,  ©  -  supplement  to  > 
cal  latitude,  u  -  coefficient  of  turbulence,  g  -  accelerat 
due  to  gravity,  E  -  gas  constant,  €  ~  flux,  of  heat  per  rna 
unit  dependent  upon  radiation  and- transition,  of  water  in  at¬ 
mosphere,  from  one  phase  state  to  the  other,  c  -  specific 
air  under  constant  pressure,  &  -  adiabatic  temperature  gra¬ 
dient. 

Boundary  conditions  and  initial  data  are  indispensible  i 
the  complete  determination  of  the  problem. 

As  boundary  conditions  with  the  earth  surface  we  take  t 
transformation  of  the  vertical  velocity  to  aero 

.  < 

w  s  0  with  z  z  0 
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and  the  conditions  of  the  free  surface  at  the  upper  limit  of 

the  atmosphere 


dp 

dt 


y  0  with  z  03  . 


(. 


Elementary  data  are  pressure  fields  and  temperatures  at  t  ~ 
p(x,  y,  s)  and  T(x?  y,  s). 

It  is  advisable  to  go  over  to  the  coordinate  system  x! , 
yt }  pt  t  ’,  where  pressure  p  is  assumed,  to  be  an  independent 
variable;  axes  x5  and  y!  are  on  the  isobaric  surfaces.  The  < 
tit  vide  of  the  isobaric  surfaces  z  is  new  considered  a  func  *.■ .<■ 
of  coordinates  x.\  yr ,  p  and  time  t*. 

The  transition  from  the  elementary  system  of  the  coord j 
■x,  y,  z,  t  to  the  new  system  z*,  y* ,  P,  t*  is  realized  thre 
the  following  changes: 
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The  result  is  that 


Pp 

3  z 

?? 

Pp  P  2 

3x 

:•  Sf5  ™  » 

5  3x» 

t;  -■ 

dy1 

= sf  K~' 

10 


In  the  new  system  the  time  derivative  of  any  element  has 
following  .form?  * 


wt. 
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(1,1 


This  shows  'that  in  the  system,  of  coordinate  x* ,  y* , 


J  ^  ^ 

value  ^48  vfe  ©lays  the  part  of  vertical  velocity;  this  vale ■ 

rl  v  * 

4,« 

gives  the  posit  loir  of  the  air  unit  in  regard  to  the  iaobaric 
surface.  The  transition  from  •  vr  to  v  is  realised  through  tr 
■correlation 
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Icuaticmis  .{1*1}  and.  1 3, *. 3 )  in  the  new  coordinate  system 'as 
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Xn  the  system  of  the  x*  »  yf ,  p»tf  coordinate  the  equate 
of  continuity  (1.4)  can  be  greatly  simplified.  To  this  end  we 
use  the  equation  of 'statics  (1.3)  &nd  correlations  (1.12)  ano 
(1.9) *  Considering  that 
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we  obtain  the  following  form  (1.4/ 
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It  can  be  easily  recognized  that  the  coordinate  system 
selected  by  us,  permits  the  use  of  the  hydrodynamic  equations 
of  the  atmosphere  in  the  same  way  as  the  equations  for  the  in 
compressible  fluid. 

How  after  the  problem  of  ascertainment  of  meteorological 
elements  was  formulated,  we  approach  the  problem  of  change  t • 
obtain  equations  suitable  for  physical  conclusions'. 


For  this  purpose  we  differentiate  equation  (1*13)  by  y* 
and  equation  (1*14)  by  x*  j  we  subtract  the  first  equation 
from  the  second.  We  assume  that  u  changes  little  in  re¬ 
gard  to  the  horizontal.  Thus  we  obtain 
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Compared  to  other  terms  of  the  equation,  components  of 

3  Re¬ 
type  ~~  •*-•*  are  generally  small  and  can  be  disregarded. 

dp  dx  * 

We  insert  the  following  designations: 
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With  the  inclusion  of  new  designations  equation  (1.20)  take 
the  following  shape: 


To  simplify  the  analysis  we  check  the  dynamics  of  atmos¬ 


pheric  motions  first 
ces.  The  latter  will 


without  the  inclusion  of  frictional  for¬ 
te  discussed  at  a  later  stage  of  the  ana' 
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For  further  simplification  of  equation  (1.22)  we  assume 
that  axis  x*  '  shows  east  and  axis  yf  north.  Then  the 
following  correlation  is  obtained: 
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« 0  -  radius  of  globe.  In  the  following,  parameter  ^  will 
always  be  considered  constant. 

Taking  into  account  all  of  the  above  conditions,  equation 
(1*22)  now  takes  the-  following  form: 
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V/e  put  down  the  initial  equations  together  with  the  limit  con¬ 
'd  it  ions  as  follows: 
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In  the  study  of  atmospheric  motion  for  which  the  eharact* 
rising  dimensions  in  regard  to  the  horizontal  are  in  th©  ord¬ 
er  1,000  kilometers,  and  the  characteristic  horizontal  veloc: 
is 

ties/10  sec/m,  the  analysis  of  equations  (1.24)  and  (1.25) 
shows  that  these  motions  may  he  considered  geostrophic  motion 

U  s  U  f  u*  , 

v  ^  V  f  v*  , 
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-  a&  component  of  gaostrophic  vrlr.ds,  u*t  v*  are  small  supple 
meriting  corrections. 

The  analysis  of  the  expression  for  the  vortex 
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shows  that  value  a  z;  i.e.  the  geostrophic  approximation  for 
the  vorthitj  includes  the  main  part  of  tiw.  VbrUoty. 

We  use  the  geostrophic  approximation  for  the  horizontal  com¬ 
ponent  of  velocity  u  and  v  in  equations  Cl. 26)  and  (1.27);  ex- 
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cluding^froa  (1.20)  with  (1.29)  we  obtain  the  following  equa¬ 
tions: 
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Together  with  limit  conditions  (1.30)  and  (1.31)  and  elemen¬ 
tary  data  this  will  suffice  to  determine  future  pressure  field; 
temperature  and  vertical  currents. 

Today  the  ascertainment  of  the  general  solution  of  the  non¬ 
linear  system  for  the  differential  equations  (1.34)  to  (1.36) 
meets  difficulties  which  cannot  be  overcome.  We  have  only  take: 
up  the  problem  to  ascertain  the  first  derivatives  ,  4-- 

and  function  z  according  to  the  given  distribution  of  the  press 
sure  fields  and  the  temperature. 

2-l  .Equation  for  Pressure  Change 

We  record  the  equation  system  (1.34)  to  (I.36)  in  the  fol¬ 
lowing  manner: 
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Index  (’)  in  the  coordinates  x,  y,  t  is  left  out. 
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and-.  P  the  mean,  pressure  on  the  surface  of 


the  earth  which  was  assumed  with  1,000  millibar. 
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We  assume  that  -parameter  m*  changes'to  a  very  small  ex~- 
tent  .  v  compared  to  the  re~;  ■ 

lative  changes  with  the  altitude  of  vertical  velocity .tr » 

From  equation  (2.2)  we  obtain  the  derivative  of  temperate. 
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V-  through  the  eouation  of  statics.  Thus  we  obtain: 
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Both  parts  of  the  equation  are  multiplied 
ferenttata  by  T .  Thus  we  obtain : 
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Excluding  derivative  — ~  from  equations  (2*1)  and  (2.5) 
vre  obtain 
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In  the  cylindrical  coordinate  system  ( r ,  Cji ,  )  where 
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altitude,  equation  (2,6)  is  changed  in  the  following  manner: 
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To  formulate  the  problem  completely,  equation  (2.£)  shouj 
be  supplemented  by  limit  conditions. 

As  the  first  limit  condition  we  assume  the  relation  direr 

ly  resulting  from  the  heat  influx  equation  (2.2}  on  the  surf, 

C&Q 

-'/el  of  1,000  millibar^on  the  basis  of  condition  (1.31)  as 
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Considering  that  T_  s  -  /  —  f — -  ,  we  put  (2.9) 
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into  the  fo.llov.ring  shape: 


r  *  A ,  \  ^ 2 

s  ^  7 


f*.  5 


•f 


=  1 


"  -  A(r ,  (j>  s  1) 


where 


A(r, f  ,  1)  a  - 


R 

g 


£c 


—  (T  .  2  )  t  -- 

l  0  f  c 


Pi 


(2, 


<X!  £ 


^  fa  “  £ 


0,1. 


g 


(2, 


The  second  limit  condition  is  obtained  from  the  heat  in¬ 
flux  equation  for  the  upper  limit  of  the  atmosphere  in  re  go¬ 
to  the  .volume  unit  of  the  air: 
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On  the  basis  of  (1.30)  v/e  have 
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Moreover,  the  quantity  of  heat  £*  flowing  at  the  volume, 
unit,  is  very  small  in  the  upper  atmospheric  layer.  Therefc;. 
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Finally,  lira  pu  *  lim  pv  »  0  and.  the  derivatives  of  tern* 
I'ifd  C*K> 

e  J&IW0&O  £// 

perature  and  ■— —  are  !"  .  .....  L  the  isobaric  surfaces 

in  regard  to  the  large-scale  motions;  as  a  result 
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In  this  way  we  obtain  the  relation 
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Condition  (2.14)  can  be  put  down  in  the  following  satis- 
ory  manner  if  we  take  (2.3)  into  account: 
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In  this  manner  the  work  done  for  the  ascertainment  of  the 
first  derivatives  according  to  time  led  from  the  altitude  of 
isobaric  surfaces  to  the  integration  of  the  inhomogeneous,  d: 
ferential  equation  of  the  second  order  (2.g)  under  limit  con* 
ditions  (2.10)  and  (2.15). 
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Now  we  proceed  to  the  problem.  (2„$),  (2. 10), {2. 15)* 
assume  that  functions  f^Cr,  f  *  £  ^  and  A(r,f  ,  1)  can  he  pre> 

sertted  in  the  following  form: 
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HAy©  JTy. (x)  is  Bessel’s  function  of  the  n-order ,  and  the 
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sponding  expressions  are  examined* 

VIo  look  for  the  solution  of  the  problem  in  the  following 
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We  subject  (2*13)  and  {2.16)  under  limit  conditions  (2.10) 
and  {2.15)  in  equation  (2.6).  We  combine  the  terms  containing 
Jn(r  )  products  with  the  same  n*  and  adapt  them  to 
zero®  We  use  correlation 
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■  As  &  result  we  obtain  the  following  equation  and  the  limii 
conditions  for  a  new  unknown  function  3n(  f  „  t ) : as  follows : 
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The  solution  of  equation  (2.19)  is  composed  of  the  general 
solution  of  the  homogeneous  equation 
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and  of  the  specific  solution  of  inhomogenous  equation  (2,19). 
The  general  solution  of  homogenous  equation.  (2.22)  has  the 
foil  aw* ng  .form: 


and  C,.,0.  -  constant  quality  depending  upon  p  . 

The  particular  solution  S*  i  P,()  of  the  inhomogenous  < 
tion  (2.19)  v?ill  be  found  through  the  method  of  variation  c: 
constants. 
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As  a  result  of  the  solution. of  system  (2,26)  wo  obtain : 
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We  integrate  (2.2?)  at  the  limits  from  zero  to  .  The  si 
Mtrarily  chosen  integration  constants  are  excluded  and  we  oh 
tain : 


.  Consequently j  the  particular  solution  of  the  inhomogenous 
equation  has.'  the  following  form: 
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In  this  manner  we  obtain  the  desired  general  solution  of  equ? 
tion  (2.19) 
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C  ,  and  C  .  are  obtained  from  the  limit  conditions  (2.20) 
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(2*21) .  The  result  of  condition  (2*21)  is  that 
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Value  C-i  can  be  determined  from  condition  (2,20) 
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On  the  basis  of  (2,30)  and  (2.31)  in  (2.29)  we  obtain 
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V/e  also  take  the  values  of  Fn(  f  t  n  )  and  G^i-y  )  from  (2.1 
as  a  basis  for  (2.33)  and  change  the  sequence Ibt  integratior. 


29 


The  result  is  as  follows: 
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If  we  take  the  vertical  line  crossing  the  zero  of  the  cc 
*.  -dinate  system  (r  s  O)1  as  a  testing  point,  we  get  n  sj4  0 
all  values  if  Ja(0)  -  0  with  n  ^  0  and  Jo(0)  -  1  arc  take 

into  aeoeunt. 

*4n)  s  0  and  M*(n)  =  0 
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-*•  We  can  always  combine  the  vertical  axis  of  the  coordinate 
the  studied  vertical* 
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Only  functions  and  M  lu'  differ  from  sera. 

In  this  manner,  the  sumtotal  { 2.34)  is  transformed  to 
composed  solution,  with  r  -  0. 
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The  result  of  (2.37)  is  that  the  pressure  change  at  the 
points  which  are  situated  on  any  vertical  line,  is  determine 
by  the  mean  value  of  the  exact  functions  of  the  pressure  a: 
temperature  fields,  in  the  circumferences  on  the  isobaric 
surfaces,  the  center  being  on  the  examined  vertical. 

This  result  can  be  obtained  through  a  simpler  mathemati¬ 
cal  process.  We  integrate  equation  (2,8)  with  f  , ,from  zero 
to  2/r  .  VJe  find  the  solution  for  the  deriv  atives  —  |  ,  ob¬ 
tained  with  the  circumference  of  radius  r  j  thus  we  assume 
that  in  the  obtained  solution  r  -  0„ 

In  the  following  this  method  will  be  applied  for  the  ob 
tainment  of  the  solution  for  and  V  . 
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To  simplify  solution  (2,37)  we  try  to  find  the  function 
of  influence  referring  directly  to  Afr’y  )«  'Por  this  purpose 
we  use  expression  -  y~  /|A(r  * ,  )  instead  of  ?  > 

(2.37).  We  integrate  in  stages  and  obtain: 
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Considering  (2.41)  and  (2.42)  we,  finally,  present  the 
solution  in  the  following  form: 

1  oo  2  rr 


{  { 


ft  2  tr 


nr 


g 


• -*  (  £  5  A  Zi)  f  ft  *"**” 

I-  '  r 


M-j_(  ^  ,  r,  r|  }rd  <pdrdyj 


o  o  o 


i  co  2?r 


4*  <■*«»>•«• 

2r 


R  (  g 

**  «w 

g 


fg  f  1 

!•“.*» (T.  z)  ,  r ,  }rdd> 

\  l  c„  j  '  ‘ 


o  o  o 


drd  . 

(2. ht 


The  influence  function  Kj  ( f  ,  r ,  yj  )  and  M^t  ?  ,  r,  ^  )  is 
characterised  "by  the  dependence  area  of  the  solution  of  the 
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meteorological  element  fields  in  the  surrounding  space.  . 

To  calculate  functions  of  K* ,  and  it  is- appropria¬ 
te’  make  certain  changes  in  (2,39),  (2.41)  and  (2.43).  The 
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Oj  1  5 

tive. 


.  t  f 


|Xj  F  are  presented  as  instructive 
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As  a  result  function  M3  ( <T  ,  r,  \  )  Is  con.erted  to  the 
algebraic  sumtotal  of  three  integrals  of  type  (4): 


x  ,  r)  m 


{  ei0»  1-  isilil  *  e  £-21^ 
)  F  V  In2  --  4  J 

o  “  X 


-  vVln2“~  f  r2 


(2,45) 


where  p  s*y  >  jp*  ,  and  x  takes  value  f>|  ,  ~'S  • 

If  we  use  the  symbols  of  (2,45)*  function  Mj{£,  r,  ) 
is  presented  in  the  following  manner: 

%(?,  *%  *|  )  *  zWr  L^1<r  (fU*  *)  *  ^  r\  * 

’ i  (2.46) 

To  calculate  M*( x  ,  r)  we  use  the  transformation 


7 + 1* 


Considering  (2.45)  we  obtain: 


r  ^  --“J™-  OX  = 


^.  {  A-  ,r)  (  -ft 

t  f--  -  I , 


X 

s* 


“  2*  <f  ( A  »  r)  dX  .  (2.47) 
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We  introduce 


<jr  (it  r) 

ms  *b>  *er.  <nt  t»f  »  ax 

V  x 


U  ( ic  ,  r ) 


(2.4 


As  a  result  we  can  present  function  (  £  ,  r,  vy  )  through 
the  following  expression: 


r  ,  r,  f]  1  »  — 

I  2 

a 


s 


;  r 

\  > 


-;r  U  /•?•»  r)  f  ™  U  fX,.  rj  f  U(C^,r) 
1  ' 


•^»r 


W  fUf  KW 


v  2 


\  4 

4  (jij  r-4  f  r_1u  <x ,  r)dr , 


(2.49 


<)  £  . 


level  for  which  **««  is  determined ,  r,h  -  the  vari 

C?  l'  t 

able  quantities  of  integration. 

To  calculate  function  N{  f  ,r,  *j  }  vre  use  the  transformat i 

Op  'S  co 

/  ']**"  qb  j  /  „  lQ(rp  ) 

v  yt.  J0*r/>  }  e  y  ‘1  j^d p  (2.50 

o  o 


and  then  also  correlations  (2.45)  and  (2.4S).  Then  N(j°  ,r, 
becomes  the  smntotal  of  three  integrals  of  type 


co 

/ 

X  (x,  r)  «  X;  2  J0(  r  p  d  p  s  / 
0 


2x 


1  Xn 

w  «*»  **w  trw  *» 

2 


r  t  ♦  Y^2  *-  #  r2) 


■  w«.rn'.iwia>wwinf>i>»iiint«ii 

p  x  <. 

In  ”  *****  r^ 

1 


U( x  .  r)  . 


(2.5: 


Taking  note  of  (2.46),  (2.48). and  (2.51)  function  M2  can 

now  be  described  in  the  following  manner: 


.  »  *  2  T  A  J  „  J.4  x  /-w  r)  -  I(Cn 

M2(C  ,  r,i|  )  s  --  |y  1  if*  J  %  ((  /  ( 


,  r)[  f 


U(f«.r)/  4 d *•)•  (2*5* 


£i  U  U  '*T 


.  We  note  particularly  that  vrith  f  s  ls  the  influence  func« 
tion  M-i  and  M2  has  a  usual  form  as  follows: 


Mx  (1,  r,  rj)  r  M*  ( *[  ,  r) 


(2 


C' 
*  J  w 


M2  (1#  r,7j  )  s  -  I  ( **}  »  2*)  +  ^  ^ 


(2*54 


Functions  U  (X,  r),  I  (*  ,  r) ,  M*  (*  ,  r)  and  the  fields 

of  functions  Mx  (f  ,  r,  ^  )  and  (C  »  >  for  <^£erent 

cutting  through  the  vertical  plane,  are  graphically  shown  in 

Fig.  1-12. 

We  will  analyse  the  reaction  of  these  functions  with  dif¬ 
ferent  values  of  the  independent  variables  X ,  r»  *f  •  ^ 

Function  U  (  x  ,  r)  (Fig.  X)  is  V  (1,  r)  =  --  £  2  with 

1,  and  u  (0,  r)  =0  with  0;  with  a  small  X  the  functi 

1  1 

tends  toward  zero  such  as  - — ■£*  *  ~v,rith  r  ~  C  U  ^  Z  *  °^ln  * 

In 
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with  growing  r  and  a  fixed  value  x  function  U( x  ,  r)  decrease 
in  a  regular  manner;  with  r  -r>  oo  it  tends  tovrard  zero,  as 

A e  "  'o'  .  The  function  shows  its  characteristic  trait  at  pc. 

r  c 

(  x «  Ij  r  k  0)* 

With  xs  1  and  r  A  0  function  I  (x  ,  r)  (Fig.  2)  equals 
zero;  withxs  0  X  (0,  r)  -  0;  with  a  small  x  the  function 


tends  toward  zero,  as  - •  V/ith  r  «  0  I  [x  ,  0)  r 

,  2  In  A 

2  «f  In  A 

— ................  With  growing  r  and  a  fixed  value  /  function' 

21n2  A 

I  (x  ,  r)  decreases  in  a  regular  manner;  with  r  oo  it 


1  -  X 

tends  toward  zero,  as  — — e  2  .  The  function  is  character- 

2r2 

istic  at  point  ( x ->  1,  r  s  0). 

The  characteristics  of  function  M*(X ,  r  )  (Fig.  3)  are 
similar  to  function  U  (X  ,  r).  Values  M*(x  ,  r)  are  1.5  to 
two  times  higher  than  values  U  (%  ,  r). 

With  x  s  0  M*(0,  r)  s  0;  with  r  -$►  oo  M*(X  ,  r)  appear? 
1  «■ 

as  —  e  (Fig.  4).  The  tendency  of  the  function  toward 


zero  is  governed  by  the  rule  -y-  . 

In  w- 

Function  (1,  r,>|  )  Fig.  5)  describes  the  integral  act 

.of  "the  dynamic  factors"  -~(z,Az)  +  y-  in  regard  to  the 

pressure  change  at  the  point  on  the  earth’s  surface.  The  act 
distance  R  =  a  750  km  corresponds  to  the  relative  length 


Fig.  I 

Graphic  Presentation 
of  Function  U(  x.  ,  r) 
with  Different  r' 
Values 


Fig*  -  2 

Graphic  Presentation 
of  Function  I{ x  ,  r) 
with  Different  r 
Values 
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%■ 


Fig.  3.  Graphic  Presentation 
of  Function  M*  (  x  ,  r)  with  Different 
r  Values 


i 

:  i  . 

j  I  .X-"/* 


U- -L - 1 J— — 1 J-~ - -  — i = - 1 - — - . - 1 

O  01  04  Ob  c>.&  iff  2,0  i  o  tor 


Fig.  4*  Graphic  Presentation 
of  Function  Mv  ( x  ,  r)  with  Varying 
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Attention  must  be  given  to  the  vast  expansion  of  the  action 
area  of  the  dynamic  factors  in  horizontal  direction.  For 
better  illustration  and  analysis  of  vanishing  of  action  on 
the  part  of  the  dynamic  factors  with  increasing  r  values, 
we  will  show  the  influence  of  rM^  (1*  r,^  },  in  a  graphical 
manner,  describing  the  importance  of  average  values  of  dyna¬ 
mic  factors  JL  (z.i>z)  t  relative  to  circumferences  of 

i  ■  '  jjt 

radius'  r  (Fig.  6).  We  see  that  function  rM^  vanishes  very 

slowly  with  increasing  r.  With  r  oo  rMj  ( 1 ,  r,>|  }  decrees 
r  ■ 

**  ;"vp 

as  £  . 


In  actual  practice,  the  influence  of  dynamic  factors  on 
the  change  of  surface  pressure  will  always  be  limited  by 
radius  R  of  the  2,000  km  order,  because  the  average  values 


g  tz 

-,**  (  2  ,  &  ?,)  f  ft  -«■- 
L  r  } t 


relative  to  the  circumferences  2,000  ~ 


3,000  km  are  actually  small. 

The  limited  effect  on  the  change  of  surface  pressure  and 
dynamic  factors  in  the  upper  layers,  can  be  seen  in  a  cleare: 
manner  in  Fig.  6.  With  an  equal  distribution  of  the  dynamic 
factors  according  to  altitude,  their  effect  from  the  surface 

limited  by  radius  r  *  1  (300  millibar)  is  reduced  tv/ofold 

\  : 

compared  to  the  effect  of  these  factors  from  the  900  milliba 
surface  which  is  limited  by  the  same  radius. 
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i 

') 


Function  M2  (1,  rf^  )  (Fig.  7)  gives  an  idea  of  the  effec 

of  thermal  factors  -S  (T,  z)  f  or,  to  be  more  precise, 

t  '-p 

of  the  effect  of  local  heat  flux  in  the  atmosphere  and  the 
change  of  surface  pressure.  The  negative  values  of  function 
M2  (1,  r, y  }  in  the  surrounding  area  of  the  atmosphere  show 
that  at  any  level  of  the  atmosphere  in  this  area,  the  local 

^hito  fo'i 

heat  flux  causes  a,  pressure  drop  at  the  earthy  heat  emission 

increases. 


Fig.  5.  Field  of  Function  Fig.  6.  Field  of  Function 

M1  (  r  ,  r,  37  )  with  1.0  rMx  (i,  r,|  ) 

Compared  to  the  action  area  of  dynamic  factors,  the  regi- 
of  activity  of  thermal  factors  appears  to  be  more  limited. 
Indeed,  the  maximum  value  of  function  M2  (1,  rt’ij  )  along  the 
vertical  line  with  r  =  1  is  more  Phan  twenty  times  smaller 
than  the  maximum  value  of  (1,  r,^  )  along  the  verical  li* 
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with  r  s  0,2,  during  the  same  period  in  which,  for  instance, 
the  function  of  influence  Mq_  {1,  r, /|  }  decreases  4  times. 

The  difference  between  the  areas  Of  influence  of  thermic 
and  dynamic  factors  can  be  recognized  by  the asymptotic  re¬ 


action  of  functions  Mg  (1,  r,  ^  }  and  (1,  r,  ). 

Function  1  ( jj  ,  r)  which  appears  as  the  main  part  of  function 
M2  (1,  r , <[  ),  tends  toward  zero,  with  higher  r  values,  as 


but  (1,  r,^j  }  shows  the  same  reaction  as 


We  also  want  to  show  function  r  Mg  {1}  r,  )  which  represent 
the  summary  influence  of  thermal  factors  on  the  circumference 
of  radius  r,  as  well  as  the  change  of  surface  pressure.  The 
analysis  of  function  .  r  Mg  (1,  r,^  )  shows  that  the  influence 
of  thermal  factors  in  the  upper  layers  of  the  atmosphere  on 
the  change  of  surface  pressure,  is- usually  very  small.  Inasmut 
as  the  average  values  of  the  thermal  factors  CT,  z  )-*-/-  on 

.......  i  ;  C  P 

circumference  r  with  a  sufficiently  high  r,  can  be  considered 

important  only  in  the  case  of  greater  thermobarometric  dis¬ 
turbances  in  the  atmosphere*  On  the  basis  of  the  graph  shown 
in  Fig.  S  it  can  be  maintained  that  the  relatively  important 
influence  of  thermic  factors  of  the  upper  layers  on  the  chang 
of  surface  pressure  can  exist  only  in  the  case  of  great  therm 
barometric  disturbances. 


43 


r 


Fig.  7«  Field  of  Function 
’  (  f  r  ,  }  with,  f  «1.0 


Fig.  S.  Field  of  Function 
rM2  (1,  r,^  ) 


How  we  'examine  the  nature  of  the  functions  of  influence 
of  (  f  ,  r,  fj  )  and  M2  (  f  ,  r,>|  }  with  (* <£  1,  i.e.  the  effec 
of  the  above  mentioned  dynamic  and  thermic  factors  on  the 
change  of  pressure  at  a  certain  level  above  the  surface  of  th 

earth* 

Fig*  9  and  10  show  functions  (0,7,  r,  ^  )  and  %(0.5, 

:r,  m  },.  In  principle,  they  do  not  differ  from  functions  which 
are  similar  to  them,  i.e.  from  functions  for  (si,  which 
warp  studied  above.  The  difference  is  only  that  the  area  of 
maximum  values  of  the  function  differed  here  from  levels 
£*  0,7  and  0*5,  that  is  to  say  from  levels  of  examined 
points#  _ 
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Fig.  9.  Field  of  Function  ■  Fig. ,10.  Field  of  Function 

Mi  (£  ,  r,rj)  with  £  -  0.7  K  ,  r,  ^  )  with  C  r  0.5 


A  completely  different  type  is  function  (C  »  r»  'j ) 
with  ^  <  1  compared  to  M2  (1,  r,  /j  ).  Fig.  11  and  12  show 
functions  M2  (0.7,  r,/|  }  and  M2  (0.5,  r,  xj  ).  he  can  see  th«' 
the  local  flux  of  heat  in  the  upper  part  of  the  atmosphere 
relative  to  level  £  causes  a  drop  of  temperature ;  on  the 
other  hand,  the  heat  flux  in  the  lovrer  part  relative  to  lev. 
£  leads  to  a  rise  in  temperature.  The  ‘  of  heat  at 
different  levels  C,  of  the  atmosphere  causes  the  reverse 
effect. 

At  a  certain  level  the  influences  of  thermic  factor,- 
of  the  upper  and  lower  layers  on  the  pressure  change  comper. 
sate  each  other  so  that  the  pressure  change  depends  only 
upon  the  dynamic  factors.  This  level  will  be  called  "mean 
level”  of  the  atmosphere. 


Theoretically,  the  existence  of  such  a  level  was  alread 
mentioned  in  I. A.  Kibel’s  (2)  analysis.  In  the  first  place 
the  isobaric  field  at  this  level,  as  it  was  shown  by  I. A. 
Kibel,  determines  the  motion  of  barometric  and  thermal  dis¬ 
turbances  in  the  vicinity  of  the  earth.  However,  empirical! 
the  existence  of  the  steering  current  was  shown  already  at 

an  earlier  date  by  S.I.  Troitzky  (5)* 

In  each  concrete  synoptic  situation  function  M2  ( S  ,  r, 
permits  us  to  determine  the  position  of  such  a  level.  In  ge 
ral,  this  level  does  not  remain  constant  neither  relative  t 
the  vast  spaces  nor  the  time,  because  in  each  practical  ca 
its  position  will  be  dependent  upon  the  distribution  of 
heat  flux  E  •  Cp  1 1 »  z )  *  £  • 


46 


If  we  admit  that  the  local  heat  flux  is  constant  relative 
to  the  vertical  line  along  the  entire  thickness  of  the  atmo¬ 
sphere,  the  mean  level,  in  this  case,  would  be  somewhere  be¬ 
tween  the  surfaces  of  700  and  500  millibar* 


The  position  of  the  mean  level  remains  dependent  upon  the 
values  of  thermic  and  barometric  distribution  in  the  atmo¬ 
sphere.  With  increasing  r  values  the  negative  values  of  func¬ 
tion  rMg  {£  ,  r,  /j  )  in  area  n  disappear  at  a  slower  rate 
than  the  positive  values  in  area  C  »  the  conclusion  is 
that  with  the  same  relative  distribution  of  the  heat  flux 
E  {  r,  )  along  the  vertical  line,  the  mean  level  will  be 
higher  in  the  case  of  great  disturbances  than  in  the  case  of 

v  t 

small  disturbances. 

With  £  -•*>  0  solution  (2.44)  for 
following  expression: 

CD  2  ff~ 


is  converted  to  the 

d  U 


<?t 


1  (  f  rF  >z‘ 

-  i  \  m2  — (z,A  z)  *  ^ -•*- 
-  o  2r  )  )  L*  T  ^ 


ua,r) 

- — — .  rd  f>  dr 


co  2 


where 


—  j~~  {T,  z)  f 

L  °P4 


U‘(l,  r) 

— — - rdf  dr,  (2. 

c  =0  4 


/ 


u*(l,  r)  a  u(i,r)  -  21(1,  r), 


(2*5- 


/  U  (1,  r)  with  r  >  0 
U*(l,  r)  - J 

|  0  with  r  a  0. 

The  functions  of  influence  ( £  ,  r,/ j  )  and  M£( f  ,  r,^  ) 

characterize  the  region  of  dependence  of  the  studied  value 

■4--  at  the  respective  point  upon  the  environmental  setting 
it 

■of  the  fields  of  meteorological  elements;  these  functions  ca; 
also  be  interpreted  on  the  basis  of  the  principle  of  revers- 
ability  In  the  following  manner. 

'  At  a,  certain  point  of  space  (X,  y,  ( )  we  place  a  single 

•"source  of  substance* * 


t  £ 


B  «  :• 


z 

(Z,AZ)  f#-— 
r  <U1 


or  A  -  -- 


g  r 

-*•*  (T  |  2. )  4^6* 

1  CP 


Accordingly,  function  (/j  ,  r,  f  )  or  M2  ( >|  ,  r,  {  ) 


•s, i,l 


fu: 


.sea  the  change  of  pressure  in  the  surrounding  space  throu/ 
influence  of  this  source.  The  positive  source  of  the  dy- 

T 

ic  substance  B  at  the  level  .'results  in  an  increase  of  pre* 
e  in  the  surrounding  area  at  all  levels,  reflected,  by  spa* 
Ctrl  on.  M3  { ^  ,  r,<f  ).  It  can  be  easily  seen  that  function 
'(  .;•)  ,  r,  £  )  is  identical  with  (£  »  r,  ). 

The  single  "thermal  source  A"  at  level  C  results  in  de- 


dreaeed  pressure .in  the  upper  part  of  the  atmosphere  accord¬ 
ing  to  function  Mg  ( n  ,  r, f  )  as  far  as  level  £  is  c oncer: 


as  well  as  a  rise  of  pressure  in  the  lower  part  of  the  at¬ 
mosphere.  Negative  sources  of  the  above  mentioned  dynamic 
and  thermic  substances  cause,  in  an  analogous  manner,  a  re¬ 
verse  tendency.  , 

Such  an  interpretation  of  the  influence  function  gives 


us  a  clear  idea  of  the  radius  of  influence  of  thermic -barome 
trie  disturbances  in  the  atmosphere  on  the  change  of  pros  sur¬ 


charges  in  neighboring  regions. 

An  analogous  interpretation  will  also  be  given  by  the  so 

lutions  for  ---  and  <¥  ,  which  will  be  described  in  the  fol- 
3  t 

lowing  paragraphs. 


3.  Equations  for  Temperature  Changes. 

The  equations  for  temperature  changes  at  different  levels 
will  be  obtained  through  the  separation  from  the  equation 
system  (2.1)  -  (2,3)  of  the  vertical  velocity  -rand  deriva¬ 
tive  . 

r)  t 

To  this  end  we  differentiate  (2.1)  with  f  .  We  obtain: 


9 

T{ 


f  ^A 

A  — — — 

.?t  / 


+ 


x  Id 


U,A  z)  +  ll  — 


dX 


The  equation  of  statist  (2.3)  results  as  follows: 

d  (  9  z\  R  1  .^T 

— — —  — — —  1  —  — ■  — -  — —  a  — — —  . 

\  r)tj  g  r 


2  2 
LJl^ 

?IT{2 


(3.1 

(3*- 
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«■*«•*  m 

K 


g  <>Z 

7  (z- Az)  *  67 


R  1 


g 


C 


7*  (T,  £z)  +  "•  ( z j  £  T) 
t  e 


,  >* 

jm  R  mm  —  m* 

•  p  a* 


(3.3 


We  incorporate  (3*2)  and  (3.3)  into  (3.1)  obtain: 


C 


g 


♦  lu*AT)  *  ^7] :  sr/?  • 

(3.4 


/r 


To  exclude  ---*  from  (3.4)  we  dissolve  the  equation  in 

d$z 

regard  to  t'  and  differentiate  same  twice  with  £  assuming, 
as  previously,  that  parameter  m2  changes  according  to  alti¬ 
tude.  We  obtain: 


r  3  x2  A  3T 

c5?*23fJ  3" 


$  *  ^ 


J-f-lT.  .)  i 
L  0 


p 


m2 12  /r 

m'W' 


(3.5 


Now  we  multiply  (3.4)  with  m2  and  subtract  from  (3.5) 
The  result  is  as  follows: 


2. 52  2.  i  ny  ---  :  t2  (X  ,  y>C  ). 


(3.f 
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f  ( x  y,  C  )  -  — —  f  ^ 
X2  v  *  y»  *  ~ 


S  .  .  £, 

--  (T,  z )  f  -~ 

^  cp 


-~  ( TtAz }  <f  -£*  {2, AT)  + 


(3.7) 


In  regard  to  the  nature  of  limit  conditions  for  equation 
(3*6)  we  assume  the  following  (see  conditions  (2.9)  and  (2.14 


9T 

mmtm 

}T 

Tt? 


c sl 


a  Q0  (x ,  y)». 


=  0, 


(3.8) 


Q0( x »  y)  =  t*  (T0»  zo}  +  “  • 

cp 


(3.9) 


The  term  (&  -  Y)  ~2  contained  in  Q0  (x,  y)  is  disregar 

5  t 

1 

ed  in  this  analysis  because  this  term  is  of  a  higher  order  . 

In  view  of  the  fact  that  equation  (3.6)  has  a  structure 
which  is  analogous  to  the  structure  of  equation  (2.6)  for  -j 


1  The  exclusion  of  term  {)£,  - Y )  contained  in  Q(x  ,  y), 

V  Z 

is  not  of  basic  nature.  It  could  also  be  retained  as  can 
be  determined  according  to  formula  (2.44)  in  the  previous  pa 
ragraph. 
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<r*lN 


m  _ 

the  solution  for  --r  as  well  as  for  ******  will  be  determined 

d  z  &  t 

by  the  mean  value  of  exact  functions  on  the  peripheries  with 
the  center  being  at  the  examined  point.  Therefore ,  we  apply 
a  simpler  mathematical  method  for  the  solution  of  equation 
(3.6) .  The  possibility  of  applying  this  method  was  already 
mentioned  above. 

We  note  equation  (3.6)  in  the  cylindrical  coordinate  sy¬ 
stem  (r,$  ,  {  )*  r  =  f  -  Polar  angle,  ^  - 

mentioned  altitude.  We  integrate  this  equation  with  f  from 
zero  to  2 ir . 

The  result  is  as  follows: 


(i.  f  2  2.  +  -1.  1.  r  -1  :  f2(r,f ),  (3.10 

K  r  dr  dr)  £  t 


2r  2j<T" 

j  ij»  n  /  1/ 

d?  -  f2(r-P*JP  | 

l  (3.1 


d  t  2  fT  )  d  t 
o 


These  are  values  —  and  f,  ascertained  on  the  circumference 

£t  4 


of  radius  r  . 

The  limit  conditions  (3.6)  are  reflected  in  the  following 
manner: 
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2T| 

}  1 1  £  -1 


C**; 


=  Q0^it 


-  0, 


at!  |=o  / 


(3.1" 


The  solution  for  -  will  be  tried  in  the  form  of  the  Fourie* 

at 


Bessel  integral 


a? 

It 


oo 


3(f  »-f  W0tr/!  fdf 


(3.1 


V7e  assume  that  function  f  2  (r ,  f  )  ar‘d  Q0(r)  can  be  present, 
ed  in  the  form  of  the  Fourie-Bessel  integral: 


oo 

r 


00 

( 


f2(r \  f2  (r',4’  )J0(pr>.)  r'dr' 


(3.1' 


CO 

r 


00 


Q0  (r)  .  Jc(Pp)f  df  Q0(r»)J0(f  r*)r*dr* 


We  take  (3.13)  and  (3.14)  as  a  basis  for  equation  (3.10)  and 
limit  conditions  (3.12)  and  use  the  known  correlation 
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ViJ 


s 


-  /J0tf  r), 


In  this  manner  we  obtain  the  following  equation  and  limit  co 
ditions  for  S(  jj  ,  £  ) : 


2.  c2  it .  » 

H  {  3C  f 


</»  .r ). 


(3.1 


(3.1 


f2  *  (  f“2(r,»C)  JD(^  r’J^dr* 


G2(  p)  s 


co 

f 


Q0(rf)J0(^  r* )  r'dr  ’ . 


y 


(3.1 


In  view  of  the  fact  that  equation  (3.15),  as  far  as  its 
structure  is  concerned,  is  analogous  to  equation  (2.19)  exa¬ 
mined  in  the  last  paragraph,  the  general  solution  of  equatio 
(3.15)  will  have  a  form  which  is  analogous  to  (2.29). 


S(  f ,  c  > 


OiC1’  t 


Now  we  assume  that  in  (3.20)  r  s  0,  As  a  result  we  obtain 

i  *p 

the  solution  for  ---  at  the  points  which  are  located  at  varic 

V  t 

levels  £  ,  along  the  coordinate  axis, 
oo 

j  f2(rl.]  >m3(C 

o  !  ^  .  .  o 

I(f  ,  r')r*dr»,  (3.2 


oo 


o 


These  are  functions  which  had  already  been  used  previously 

We  take  f2(rS?|  )  from  (3.11)  and  (3.7)  as  a  basis  for 

(3.21)  and  partially  integrate  thermic  factor  &(T,  z)  /  --- 

L  cp 

Thus  we  can  finally  give  a  description  of  solution  (3.21)  for 

•4*%  as  follows  1  :  ,  ,  ,  ;; 

?z 

1  Derivative  can  be  directly  obtained  from  solution  (2.3 
for  y|  by  differentiating  with  C  and  utilizing  correlation 


oo 


1 


)r1drt&^  + 


Q0(r* ) 
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ca  2jr 

3?  Iff 

c?t  ’  2r  )  J 

o  o 


2  (T  ,  z  )  f  -~2 

l  o>  o  7  Cj> 


I  { f ,  r )  rd  f  dr  4- 


1  cd  2r 

,  -i-  (  (  (  ,  ’ 


2rf 


2r> 


o  o  o 

1  co  2  rr 


~ (T ,  2 )  f  ---  j  M^(  f  ,r ,  )  rd  <p  drd  ^  4 


~7 


2?r 


r 


m 


"  ~  ( T »  3 }  f  -j*  ( z ,  A  T )  f  ^ 


c>X 


,  *v? 


drd  a  , 


I 


(3.2- 


M4  (f  »  r>)l  )  = 


>  M, 

X[“ 


(continued) 


Then  we  obtain: 

oo  2/r 

i  7 


2! . .  f.  CjL  Lti 

7 1  "  R  ^  31’  D  ty 


2t  ”  2r 


QQ(r,tp )Gi.(f >r)rdfdr  4 


o  o 


O  O  0 


gm 

R 


g  <?z 

“£*  (  Z  ,  A  Z  )  -f 


G2(f,r^)rd^drd  ^ 

.  ) 


a,  7 


*  ~  ( T ,  z )  4  - 
*■  c 


p 
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The  functions  of  influence  and  are  expressed  by 
functions  of  type  U(*  *  r)  and  li*  ,  r)  which  were  used  in 
the  preceding  paragraph. 


=  7  If 


L i  o 

If 

%r 


1 

2 


M3(f 


Fig.  13-17  show  the  function  of  I(£,r)  an^  bhe  sections 
in  the  vertical  field* of  functions  M4(f,r,?j)  and  M-j(f,r,|). 

The  analysis  of  the  values  of  function  ll£>»r)  (Fig.  13) 
shows  that  the  first  term  on  the  right  (3*24)  calculated  ac 
cording  to  altitude,  steadily  vanishes  with  the  vaning  f  . 


(continued) 

G2(f»r,^)  s  - 

-1  * 

+  --  M3  {$, r,?j)  f«<M  (£/j,r), 

G1(^,r)  sf  =  K5, r)  -«*M*(?,r). 

»  C 
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Function  M^(.f,r,^)  (Fig.  14  and  15)  is  used  for  calculat¬ 
ing  the  irregular  distribution  of  the  local  heat  flux  along 
the  vertical.  The  first  two  integrals  on  the  right  (3.24) 
give  complete  data  on  the  influence  of  local  heat  flux  distri¬ 
buted  over  the  entire  atmosphere  and  on  the  temperature  change 
at  a  fixed  point.  - 

Analytically,  function  M^j^r,^)  consists  of  two  parts: 

1 

M4(f»r»^  5  x  (£>r,^)  -  —  M3  (£,r,^),  (3.26 


-Wf 


-  'I  1  (-- .  r)  -  I(f)j.r) (,  .  (3.2 


7  U 


Fig.  13.  Graph  of  Function  I(r,r) 
for  various  £  . 


In  the  environment  of  radius  r«r  0.5  X(/,r,/j)  forms  the 
main  part  of  values  M^(£,r,^).  On  the  other  hand,  A.  M^(f  ,r.* 
appears  to.  be  relatively  small.  With  higher  values  of  r{r^£l', 


the  corresponding  equilibrium  of  both  components  contained  in 
M4(r,r,«),  is  outbalanced;  but  each  of  these  components  is 


small. 

Function  ,r ,fj)  (Fig.  16  and  17)  gives  an  idea  of  the 

sphere  of  influence  of  the  dynamic  factors  f 

g  1  Li 

f  *~U,AT)  *  p  TxJ  °n  c^‘anSe  temperature  at  the  re 

spectlve  point  on  level  £  .  This  function  has  its  maximum 

values  in  the  immediate  vicinity  of  the  examined  point  and 

vanishes  with  the  increasing  r  but  also  in  the  direction  of 

the  upper  and  lower  boundaries  of  the  atmosphere.  =  0  and 


^  r  1  M3 ( f ,r )  -0,  if  £  o4  0.  ,  f 

It  can  be  easily  recognized  that  with  higher  r  values 
function  M3(^,r,»)  is  presented  in  an  asymptotic  manner  as 
follows:  “2“  ^  **  ^  •  According  to  the  increase  of  r  values, 
the  area  of  maximum  values  of  function  with  level 

£  is  formed  at  the  upper  part. 

The  result  of  the  analysis  of  the  properties  of  function 
is  that  the  temperature  change  at  any  level  £  ex¬ 
pressed  by  the  last  term  in  (3*24)*  is  determined  by  the  dy¬ 


namic  factors  of  the  atmospheric  layers  bordering  level  { 
to  a  greater  extent.  In  regard  to  levels  the  dynamic  factor: 

_  tve  i  6-  AT~ 

of  the  upper  layers  have  a  higher  than  simi 

lar  factors  in  lower  layers. 


Nov/  we  assume  in  what  manner  the  solution  for 
cribes  the  change  of  temperature  at  the  upper  border 
atmosphere. 


Fig.  16.  Field  of  Function 
M3IC,  with  £  =  0.7 


•  Fig.  17.  Field  of  Function 
Mj(  C  >  r,  n  )  with  £a  0.5 


With  0  function  I(,f>  r)  0,  M^(  f  ,  r,/|  )  is  con- 
C  r  X  -i 

verted  to  1(1,  r)  -  —  0(1,  r)i  M-^(  ^  ,  r,«  )  and 

2.C  2 


.%(  f  *  r,  /J  )  and 


--1  U{  1.  r).  Then  solution  (3.24)  takes  the  following  shape 
with  X-y*  0: 

1  (  f  3  fg  >1  U*{l,r) 

s  - - \  {  11  ---  A.  - - --rdf dr  f 

<Jtj  /ao  2/r  j  )  [2*1  [L  c'p  fab  4  * 

ft  ft  »  ^ 


0  o 


1 

f  “ — 

2  r 


oo  2r 


o  o 


m2  Az)  +  --(z,£T)4  ^--- /  - —  rd  f  dr, 

L  i  i  *>XJ  g  =0  2  4 


0(1, r) 


Analogous  to  (2.56) 


UT (l,r)  s 


U(l,r.)  with  r  >  0, 
.0  with  r  s  0. 


(3.22 


Now  we  recall  correlation  (3.3) 


?j-IjU,Az)  =  -  -  -W.At)  i  \-(z„ 


Assuming  that 


lim 

£•»  O 


K' 


g  „  J)z 

T<MzIM-v 


-  1 
=  o  , 


(3.29 


^  Qualification  (3.29)  means  that  the  kinetic  energy  of  the  u 
of  mass  — -  according  to  altitude  does  not  increase  faster  1 
the  increase  of  «  . 

t  .  •  - 


we  obtain 

T 

~(T,^z)  i  -£{z,A1)  f  ^  “  “7 

\  »j>  ’  S  <j 

The  result  of  (3.30)  is  that  derivatives  y-  and  at 

the  upper  border  of  the  atmosphere  equal  zero;  this  means  tha 
the  temperature  on  the  surface  fa  const,  is  constant  with 
£-•>0.  In  other  words,  . 

<T,  *)f.0  =  °*  13,3 

We  further  assume  that  £  0 


=  0. 


=° 


(3.30 


1  V 


-  0. 


(3.3 


aO 


Meeting  conditions  (3.30)  to  (3*32)  from  (3.28)  we  obtain  U 
result:  :  ■ . ;  •  :  ;■ 


?Tj 

}  H  i  =° 


=  p. 


(3.3 


L.  Equations  for  Vertical,  Velocity, 

Vertical  motion  appears  as  a  component  of  the  mechanism 
of  atmospheric  circulation.  The  redistribution  of  kinetic,  r 
tential  and  internal  energies  of  air  from  one  level  to  the 
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other  is  realized  through  vertical  currents.  Therefore,  in  the 
study  of  atmospheric  processes  the  examination  of  vertical 
motion  should  be  given  careful  attention. 

At  the  same  time,  vertical  motions  are  of  interest  as  sue 
because  they  are  the  main  factor  in  the  process  of  formation' 
of  cloud  conditions  and  precipitation.  - 

To  obtain  an  equation  for  vertical  motion  £  in  the  at- 

'  '  '  '  i  r- 

nosphere,  it  is  indispensible  to  exclude  derivatives 

and  ---  from  system  (2*1)  and  (2.3).  ,  . , 

jt  (2.1) 

For  this  purpose  we  differentiate  the  equation^with  L 

■  v':  "  )  ■ 

■We  obtain: 

X  -%=  £  (a£)  4  fiu.A.)  ♦  *4?'  ■  U 

pg  n  ?f  (  at/  |t  ¥\ 

The  equation  of  statics  gives  us  correlations 


}  '  /  ?A 

7f  K;  = 


R  1  <)T 

— »  mm  mm  mm  mm  A  mm  mm  mm 

•  g'  r  ^t 


(4.2 


h  g  2  Ri  g  g 

- —  -.-(2,42)  4  ““  yf  ~“(^>2^2)  4  -r(2,AT)  4 

Jf  l  rV*  g  C  l  1 


We  take  (4.2)  and  (4.3)  as  afisis  for  (4.1).  ?hen  we  obtair. 


_L  XX  =  -  -1-  A— -  -  V  fv(T»A  z)  +  rU*AT)  4 


i  fe 


iT 


PR  M 


I  ^  ju 


<5X 


(4.* 


We  exclude  the  derivative^of  (4.4)  through  the  equation  of 
heat  flux  (2.2)  and  obtain: 


(4. 


S 


«  -  p 


--  /  A  A  (T,zH  — “  f  A(T,  Az)  i  y{z,AT, 

r  l  i  V  1  L 


^T 


(4. 


We  assume  that  the  changes  of  value  T($  -y)  in  regard 
the  vertical  and  horizontal  are  small  compared  to  the  rela¬ 
tive  changes  of  vertical  velocity.  Thus  equation  (4.5)  can 
formulated  in  the  following  manner:  ... 


£2  Xt-X  i  m2/ir 


m2  -  parameter  introduced  above. 

The  limit  conditions  for  equation  (4.7)  are  as  follows; 


69 


P--V— 


tT0  a  g^0  -  small  quantity. 

We  now  come  to  the  cylindrical  coordinate  system  r, f  ,  (  t 
already  used  above. 

In  this  coordinate  system  equation  (4.7)  assumes  the  fol¬ 
lowing  shape :  . .  r  :v'; 


?  1  ^  1  YT  a 

^  f  **"  *“•**  t  — -  ^  s  -  ,  f*  ,  f  )•  (4.9) 

r  dr  Jr  r 2  *■*  1 


Integrating  this  equation  with  f  from  zero  to  2 Jr  ,  we 


2/J- 


obtain  the  following  equation  for  function  v  z.  tdf  ./> 


r2-Sf  ~£'£"  CF^’V 


2r 


F3(r,f  )  « 


2r 


d*3(r, f  ,  C)  4f 


(4.10) 


(4.11) 


The  limit  conditions  of  equation  (4.10)  are  as  follows: 


— irl 


r  0» 

!gs° 


si 


=  rft(r). 


(4.12) 


The  solution  of  equation  (4.10)  is  obtained  through 

<s> 

V(r,X)  = (  Stp  )  J0(rp  (4.13 
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As  we  have  done  above,  we  assume  that  function 
and  /F0(r)  can  be  presented  as  intervals  according  to  Fourie- 


Bessel 


co  co  - 

?;£(r>C)  .*  ^  J0(rj5)  <j  dj>  ^  ^r’Jr’dr’, 


To<r)  5  \  Jo^rf^  f>&p  l  trQ(r9)  J0(f>r')r* dr*  . 


(4.14) 


V/e  take  (4.13)  and  (4.14)  as  basis  for  (4.10)  and  (4.12). 
Analogous  to  the  above  we  obtain  the  following  differential 
equations  and  limit  conditions  for  the  new  unknown  function 

S<f  .f)! 


\2S 

c  7:2  -  f2s  =  <°f3 »r>'» 


(4.15) 


s|  .0, 

iC»o 

^  ^  r  ^*3  ( y°) 


(4.16) 


Following  is  the  general  solution  of  the  equation  (4.15) 


c1?5- '/  *  -L  j  Fj^^iyr 


In  the  selection  of  and  we  will  see  to  it  that  solution 
(4. IS)  meets  limit  conditions  (4.16).  The  result  is  as  follows 


Bearing  in  mind  that  and  G^C^)  ca?l  ©^pressed 

the  integral  way  (4*17) *  we  will  use  it  as  a  basis  for  (4*19 
Moreover,  the  obtained  expression  is  taken  as  a  bases  for 
(4.13) .  Now  we  change  the  order  of  integration  and  assume  th. 

r  s  0. 


Fig,  IS,  Graph  of  Function  £l(S,  r.) 
for  various  £  , 


Then  we  obtain : 

'1  CD 


r($)  «  -  ^  j  F3(r\^)  tM3(£,  r*^)r*dr*d^  * 


o  o 


CO 


+  .tz<r9){  KC,  r«)  r»dr», 


(4.^ 
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M^(f,  r*,^)  and  I(J,  rM  -  Functions  (3*22)  and  (3*23 )  in¬ 
troduced  previously. 

Substituting  F-j(r,^)  and  f0(v)  in  (4.20)  according  to 
(4.11)  and  (4.12)  through  (4.6),  we  formulate  solution  (4.20) 
in  the  final  form  as  follows: 

1  co  2 /r 

r(  {  )  „  -L  f  |  |  P  _.j  .  i.  A  (T , x)  +  (T,A  z)  f  (z,AT) 

o  o  o 

f  6L  + 

cp 

oo  2*r 


o  o 

The  graphic  presentation  of  function  f I(  f,r)  with  variou 
£  (Fig.  IS)  shows  that  the  second  terra  in  the  right  part  (4.2 
quickly  vanishes  at  the  respective  altitude. 

Considering  that  a  g^p0  usually  does  not  exceec 
10  millibar  (12  hours),  we  may  disregard  the  second  component 
in  (4.21). 

Fig.  19  and  20  show  the  fields  of  functions  CM^(f,r,^) 
for  different  £  . 


(  Cm3(£\  r,^)rdfdrd>^  + 


74 


Fig.  19.  Field  of  Function 
CM3 (£,  r^)  with  0.9 


Fig,  20.  Field  of  Function  £,WU(£,r>w) 
with  0,7  ( 
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show  that  In  the 


The  properties  of  function  M^,  rt)j) 
formation  of  vertical  motions  in  the  intermediary  tropospher 
the  dynamic  processes  in  the  intermediary  and  upper  parts  of 
the  troposphere  play  the  decisive  part.  The  vertical  motions 
at  the  level  of  3-5  km  are  mainly  determined  by  particularity 
of  the  pressure  fields,  temperature  and  heat  flux  at  the  3-£ 
levels. 

Equation  (4.21)  shows  that  in  the  source  region  of  the 
heat  flux  the  terms  of  type  A  (?,z)  +  wil1  always 

provide  the  anabatic  motions..  But  in  the  heat  discharge  re¬ 
gions  they  will  provide  the  catabatic  motions*  Consequently 
part  of  the  heat  influx  entering  any  important  region  of  th< 
atmosphere  is  transferred,  together  with  the  vertical  curret 
to  higher  layers. 

On  the  other  hand,  the  heat  discharge  is  partly  compen¬ 
sated  by  the  transfer  of  heat  through  vertical  currents  fror 
higher  layers. 

This  points  out  the  important  role  played  by  vertical  me 
tions  of  the  atmosphere  in  the  redistribution  of  thermic  ene 

Factor  in  the  solution  for  't'  which  appears  ir 

i 

the  initial  equations  of  the  parameter  change  of  the  coriol- 
according  to  latitude,  is  important  only  in  the  presence  of 
great  thermic  disturbances  of  the  atmosphere.  In  the  easter. 
part  of  the  thermal  crest  this  factor  becomes  an  anabatic  c 
ponent  of  vertical  velocity,  in  the  western  part  of  the  ere 
however,  a  catabatic  component. 
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